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A Maple package for computing Grobner bases of linear difference ideals is described. The underlying algorithm 
is based on Janet and Janet-like monomial divisions associated with finite difference operators. The package can 
be used, for example, for automatic generation of difference schemes for linear partial differential equations and 
for reduction of multiloop Feynman integrals. These two possible applications are illustrated by simple examples 
of the Laplace equation and a one-loop scalar integral of propagator type. 
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1. INTRODUCTION 

As it is argued in 1 , Grobner bases form the 
most universal algorithmic tool for reduction of 
loop integrals determined by linear recurrence 
relations derived from the integration by parts 
method. These recurrence relations can be con- 
sidered as generators of an ideal in the ring of fi- 
nite difference polynomials. Hence, to determine 
the minimal set of basis (master) integrals and to 
express other integrals in terms of basis integrals 
it suffices to compute any Grobner basis, i.e., for 
any ranking (ordering) on difference operators. 

Another important application of difference 
Grobner bases is generation of difference schemes 
for partial differential equations |2I3| . A dif- 
ference scheme is constructed by elimination of 
partial derivatives from the system of difference 
equations obtained by discretization of the dif- 
ferential equations and certain integral relations. 
Therefore, the corresponding elimination ranking 
is to be used. 

In this paper we present a Maple package for 
computing Grobner bases of linear difference ide- 
als. The package is a modified version of our 
earlier package |4] oriented towards commuta- 
tive and linear differential algebra and based on 
the involutive basis algorithm p]. The modified 
version is specialized to linear difference ideals 



and uses both Janet and Janet-like divisions |S] 
adopted to linear difference polynomials |7|. 

2. THE MAPLE PACKAGE LDA 

2.1. Outline of the package 

The package LDA (abbreviation for Linear 
Difference Algebra) implements the involutive 
basis algorithm 3 for linear systems of differ- 
ence equations using Janet division. In addition, 
the package implements a modification of the al- 
gorithm oriented towards Janet-like division |H] 
and, thus, computes Janet-like Grobner bases of 
linear difference ideals jjj. 

Tabled collects the most important commands 
of LDA. Its main procedure JanetBasis converts 
a given set of difference polynomials into its Janet 
basis or Janet-like Grobner basis form. More pre- 
cisely, let K be the difference ring [Jj of polynomi- 
als in the variables 6^oy k , /i £ Z> , k = 1, . . . , m, 
with coefficients in a difference field K containing 
(Q for which the field operations can be carried 
out constructively in Maple. We denote by M.l 
the set of linear polynomials in M. Given a gen- 
erating set F C Ml for a linear difference ideal / 
in R, JanetBasis computes the minimal Janet(- 
like Grobner) basis J of I w.r.t. a certain mono- 
mial order (ranking). The input for JanetBasis 
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3 However, Janet trees for the fast search of involutive divi- 
sors are not yet implemented since Maple does not provide 
efficient data structures for trees. 
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consists of the left hand sides of a linear system 
of difference equations in the dependent variables 
y , . . . , y m , e.g. functions of x-y, . . . , x n . The 
difference ring R is specified by the lists of inde- 
pendent variables Xi, ■ ■ ■ , x n and dependent vari- 
ables given to JanetBasis. The output is a list 
containing the Janet (-like Grobner) basis J and 
the lists of independent and dependent variables. 

After J is computed, the involutive/JT— normal 
form of any element of Ml modulo J can be 
computed using InvReduce. Given p € Ml 
representing a residue class p of the difference 
residue class ring M./I, InvReduce returns the 
unique representative q 6 Ml of p which is not 
involutively/j7— reducible modulo J. A K-basis 
of the vector space Ml/(/ fl Ml) is returned by 
ResidueClassBasis as a list if it is finite or is 
enumerated by a formal power series |1 0| in case 
it is infinite. For examples of how to apply these 
two commands, cf. Section 

Given an affine (i.e. inhomogeneous) linear sys- 
tem of difference equations, a call of CompCond 
after the application of JanetBasis returns a 
generating set of compatibility conditions for the 
affine part of the system, i.e. necessary conditions 
for the right hand sides of the inhomogeneous sys- 
tem for solvability. 

Moreover, combinatorial devices to compute 
the Hilbert series and polynomial and function 
etc. H] are included in LDA. 

For the application of LDA to the reduction 
of Feynman integrals, a couple of special com- 
mands were implemented to impose further re- 
lations on the master integrals: By means of 
AddRelation an infinite sequence of master inte- 
grals parametrized by indeterminates which are 
not contained in the list of independent variables 
is set to zero. Subsequent calls of InvReduce and 
ResidueClassBasis take these additional rela- 
tions into account (cf. Section 0J). 

Finally, LDA provides several tools for deal- 
ing with difference operators. Difference opera- 
tors represented by polynomials can be applied 
to (lists of) expressions containing y 1 , y m 
as functions of the independent variables. Con- 
versely, the difference operators can be extracted 
from systems of difference equations. Leading 
terms of difference equations can be selected. 



2.2. Implementation details 

We consider difference rings containing shift 
operators which act in one direction only. If a 
linear system of difference equations is given con- 
taining functions shifted in both directions, then 
the system needs to be shifted by the maximal 
negative shift in order to obtain a difference sys- 
tem with shifts in one direction only. However, 
LDA allows to change the shift direction globally. 

Unnecessary computations of involutive reduc- 
tions to zero are avoided using the four involutive 
criteria described in [H], [HI- Fine-tuning is possi- 
ble by selecting the criteria individually. 

The implemented monomial orders/rankings 
are the (block) degree-reverse-lexicographical and 
the lexicographical one. In the case of more than 
one dependent variable, priority of comparison 
can be either given to the difference operators 
("term over position") or to the dependent vari- 
ables ( "position over term" / elimination ranking) . 

The ranking is controlled via options given to 
each command separately. The other options de- 
scribed above can be set for the entire LDA ses- 
sion using the command LDAOptions which also 
allows to select Janet or Janet-like division. 

3. GENERATION OF FINITE DIFFER- 
ENCE SCHEMES FOR PDES 

We consider the Laplace equation u xx +u yy = 
and rewrite it as the conservation law 

<j) —Uydx + u x dy = 0. 
Adding the integral relations 




u x dx = u(x j+2 ,y) -u(xj,y), 



rVk+2 

I Uydx = u(x, y fe+2 ) - u(x, y k ) 

J Vk 

and using the midpoint integration method we 
obtain the following discrete system: 

( -{8 x -8 x 8 2 y )ou y + (8 2 x 8 y -8 y )ou x = 0, 

I 2Ah6 x ou x - (61 -l)o U = 0, (1) 

[ 2Ah6yOu y -(6l-l)au = Q, 

where 8 X and 8 y are the right-shift operators 
w.r.t. x and y, e.g., 8 X o u y (x,y) = u y (x + l,y). 
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Table 1 

Main commands of LDA 



JanetBasis 

InvReduce 

CompCond 

HilbertSeries etc. 
Pol2Shift, Shift2Pol 



Compute Janet (-like Grobner) basis 

Involutive / J— reduction modulo Janet (-like Grobner) basis 
Return compatibility conditions for inhomogeneous system 
Combinatorial devices 

Conversion between shift operators and equations 



Some interpretations of commands for the reduction of Feynman integrals: 



ResidueClassBasis 

AddRelation 

ResidueClassRelations 



Enumeration of the master integrals 

Definition of additional relations for master integrals 

Return the relations defined for master integrals 



We show how to use LDA to find a finite dif- 
ference scheme for the Laplace equation: 

> with(LDA) : 

We enter the independent and the dependent vari- 
ables for the problem (ux > uy > u): 

> ivar:=[x,y]: dvar : = [ux ,uy ,u] : 

Next, we translate Q into the input format of 
JanetBasis. Note that one can in general use 
AppShif tOp to apply a difference operator given 
as a polynomial similar to the ones in Q to a 
difference polynomial. 

> L: = [2*h*ux(x+l,y)-u(x+2,y)+u(x,y) , 

> 2*h*uy(x,y+l)-u(x,y+2)+u(x,y) , 

> 2*h*(ux(x+2,y+l)-ux(x,y+l))+ 

> 2*h*(uy(x+l,y+2)-uy(x+l,y))] : 

Then we compute the minimal Janet basis of the 
linear difference ideal generated by L w.r.t. a 
ranking which compares the dependent variables 
prior to the corresponding difference monomials 
("position over term" order; this ranking is cho- 
sen when using the option 2 as below). The least 
element of this Janet basis is by construction a 
difference polynomial which does not contain any 
monomial in ux and uy because ux > uy > u. 

> JanetBasis (L,ivar, dvar, 2) [1] [1] ; 

u(x + 4, y + 2) - 4u(x + 2, y + 2) + u{x, y + 2) 
+u(x + 2, y + 4) + u(x + 2, y) 

The computation takes less than one second of 
time on a Pentium III (1 GHz). 



Dividing this difference polynomial by 4/i 2 we 
obtain the following finite difference scheme: 

D 2 j (u jk )+D 2 k (u jk ) = 0, 

where 

n 2, s Uj+2 k - 2Uj k + Uj-2 k 

D 3 M = ^2 

and 

n 2, v Uj fc+2 ~ 2Uj k + Uj k-2 

are discrete approximations of the second order 
partial derivatives occurring in Laplace's equa- 
tion. 

4. REDUCTION OF FEYNMAN INTE- 
GRALS 

In order to demonstrate how to use LDA for 
the reduction of Feynman integrals, we consider 
a simple one-loop propagator type scalar integral 
with one massive and another massless particle: 

ffl, \ - T J— f 

j[K,n) .— if. n — . , / , . 

(Here k, n are the exponents of the propagators.) 

The basis integrals for this example and the 
corresponding reduction formulae were found and 
studied by several authors (see, e.g., |11I12| L 
Here we apply the Grobner basis method, as im- 
plemented in LDA, directly to the recurrence re- 
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lations which have the form: 
[d-2k - n-2m 2 kl+- 

n2+(l- -q 2 + m 2 )] f(k + 1, n + 1) = 0, . . 
[n-k-kl + (q 2 + m 2 -2-)- [> 
n2+(l- - q 2 + m 2 )] f(k + l,n + 1) = 0, 

where 

^/(fc.n) = /(fc±l,n), 2±/(fc,n) = /(fc,n±l). 
In addition, it is known that 

f(k + i,n + j)=0 Vi<0 Vj (3) 
which we will take into account later. 

> ivar: = [k,n]: dvar: = [f]: 
We enter the recurrence relations (J2J: 

> L : = [ (d-2*k-n) *f (k+1 ,n+l) -2*m~ 2*k* 

> f (k+2,n+l)-n*f (k,n+2)-n*(m~2-q~2)* 

> f(k+l,n+2), (n-k)*f (k+l,n+l)- 

> k* (q"2+m"2) *f (k+2 ,n+l) +k*f (k+2 ,n) - 

> n*f (k,n+2)-n*(m~2-q~2)*f(k+l,n+2)] : 

> JanetBasis (L , ivar , dvar ) : 

Again, the computation time is less than one sec- 
ond. Now, the master integrals are given by: 

> ResidueClassBasis (ivar , dvar) ; 

[f(fc, n), f(fc, ra + 1), f(fc + l, n), 

f(fc, n + 2), f(fc + 1, n + 1), f(jfe + 2, n)} 

© implies additional relations on the master in- 
tegrals. (Here, j is recognized as not being con- 
tained in ivar and thus serves as a parameter to 
define the additional relations.) 

> AddRelation(f (k,n+j )=0 , ivar , dvar) : 

The list of master integrals now becomes: 

> ResidueClassBasis (ivar , dvar) ; 

[f(jfc + 1, n), f(fc + 1, n + 1), f(fc + 2, n)] 

Next, we recompute the Janet basis for m = 0: 

> m:=0: J:=JanetBasis(L,ivar,dvar) : 

For the special case where m = 0, we impose the 
relation f(k + i, n) = for all i: 



> AddRelation(f (k+i ,n)=0 , ivar , dvar) : 

Now, we are left with one master integral: 

> ResidueClassBasis (ivar , dvar) ; 

[f(fe + l, n + 1)] 

We reduce f(k + 2, n + 3) modulo J taking also 
the additionally imposed relations on the master 
integrals into account. (Here, the option "F" lets 
InvReduce return the result in factorized form.) 

> InvReduce(f (k+2,n+3) , J, "F") ; 

-((2n + 4- d + 2fc)(2n + 2- rf + 2fc) 
(2 k + n - d) (n + 3 - d + k) 
(n + 2 - d + k) i(k + 1, n + l))/((n + 1) 
(2 n - d + 4) n q 6 k (d - 2 k - 2)) 

Using ResidueClassRelations one can display 
the relations imposed on the master integrals: 

> ResidueClassRelations (ivar, dvar, [i,j]); 

[f(fc, 7i + j), £(k + i, »)] 

The difference operators occurring in the last re- 
sult can be extracted as polynomials in 5k, S n : 

> Shift2Pol('/., ivar, dvar, 

> [delta [k] , delta [n]] ) ; 

[S n j , 41 

5. CONCLUSION 

We presented the Maple package LDA imple- 
menting the Janet/ Janet-like division algorithm 
for computing Grobner bases of linear difference 
ideals. Generation of difference schemes for linear 
PDEs and reduction of loop Feynman integrals 
are important applications of the package. 

These two kinds of applications were illustrated 
by rather simple examples. The first difference 
system (discrete Laplace equation and integral re- 
lations) contains two independent variables (x, y) 
and three dependent variables (it, u x ,u y ). The 
second system (recurrence relations for one-loop 
Feynman integral) also contains two independent 
variables/indices (k,n), but the only dependent 
variable /. The second system, however, is com- 
putationally slightly harder than the first one be- 
cause of explicit dependence of the recurrence 
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relations on the indices and three parameters 
(d, m 2 , q 2 ) involved in the dependence on indices. 

Dependence on index variables and parame- 
ters is an attribute of recurrence relations for 
Feynman integrals. Similar dependence may oc- 
cur in the generation of difference schemes for 
PDEs with variable coefficients containing pa- 
rameters. Theoretically established exponential 
and superexponential (depends on the ideal and 
ordering) complexity of constructing polynomial 
Grobner bases implies that construction of differ- 
ence Grobner bases is at least exponentially hard 
in the number of independent variables (indices). 
Besides, in the presence of parameters the volume 
of computation grows very rapidly as the number 
of parameters increases. 

Thus, for successful application of the Grobner 
basis technique to multiloop Feynman integrals 
with masses and to multidimensional PDEs with 
multiparametric variable coefficients we are go- 
ing not only to improve our Maple code but also 
to implement the algorithms for computing Janet 
and/or Janet-like difference bases in C++ as a 
special module of the open source software avail- 
able on the Web page http : //invo . j inr . ru 
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